Transport signatures of correlated disorder in a two-dimensional electron gas 



T. HeinzelQ, R. Jaggi, E. RibeiroQ, M.v. Waldkirch, and K. Ensslin 
Solid State Physics Laboratory, ETH Zurich, 8093 Zurich, Switzerland 

S. E. Ulloa 

Department of Physics and Astronomy, Ohio University, Athens, OH 45701-2979, USA 

G. Medeiros-Ribeircfl, and P. M. Petroff 

Materials Department, University of California, Santa Barbara, CA 93106, USA 

(February 1, 2008) 

We report electronic transport measurements on two-dimensional electron gases in a Ga[Al]As 
heterostructure with an embedded layer of InAs self-assembled quantum dots. At high InAs dot 
densities, pronounced Altshuler-Aronov-Spivak magnetoresistance oscillations are observed, which 
indicate short-range ordering of the potential landscape formed by the charged dots and the strain 
fields. The presence of these oscillations coincides with the observation of a metal-insulator transi- 
tion, and a maximum in the electron mobility as a function of the electron density. Within a model 
based on correlated disorder, we establish a relation between these effects. 



PACS numbers: 73.40.-c, 73. 21. La, 72.15.Rn 

The optical and electronic properties of self-assembled 
quantum dots (SAQDs) have been widely studied in the 
past few yearsJlTQ These experiments have been per- 
formed on individual SAQDs or on two-dimensional dot 
arrays. Recently, however, substantial progress in fabri- 
cation technology hasJed to self-ordered gcowth of such 
quantum dots in twoHB and even threa!3 dimensions. 
The ordering has been detected by scanning force mi- 
croscopy as well as by X-ray diffraction. These stud- 
ies suggest that designing electronic minibands by self- 
ordered growth of SAQDs may be feasible. 

Here, we report experimental evidence in electronic 
transport, that InAs SAQDs embedded-at the heteroin- 
terface of a Ga[Al]As heterostructureJlJ are capable of 
creating a potential landscape with substantial spatial 
and amplitude correlations. In our samples, the SAQDs 
have been grown on an atomically flat GaAs [100] surface. 
The charged SAQDs, combined with the induced strain 
fields, act as repulsive scatterers in a two-dimensional 
electron gas (2DEG). At sufficiently high dot densities, 
several effects are observed to occur in concert: (1) Mag- 
netoresistivity measurements show pronounced, periodic 
oscillations, which are interpreted as Altshuler-Aronov- 
Spivak oscillations. In addition, under certain parametric 
conditions, a resistivity minimum around zero magnetic 
field is observed, superimposed with the weak localiza- 
tion peak. (2) We observe a maximum in the electron 
mobility as a function of the electron density. (3) The 
temperature dependence of p xx as a function of electron 
density shows a icansition from metallic behavior to insu- 
lating behavior,!! 2 ] an effect that has been interpreted as 
a signature of a metal-insulator transition in two dimen- 
sions. We emphasize that in the samples under study, 
these transport features always occur together. As we 
describe below, this complex experimental behavior is a 
signature of the potential correlations^ in the 2DEG pro- 



vided by the SAQD potential landscape. 
The samples have been grown by molecular beam epi- 
taxy. A layer of InAs SAQDs has been embedded in the 
GaAs buffer layer 3 nm below a GaAs — Al0.3Ga0.rAs - 
heterointerface, at which the 2DEG resides. The density 
riQD of the InAs SAQDs varies across the wafer. ?iq£> has 
been determined, within a factor of two, by transmission 
electron microscope (TEM) studies, as well as by low- 
temperature transport experiments. For details of the 
sample growth and characterization, see Ref. |ll|. The 
wafer was cleaved into pieces, and conventional Hall bar 
geometries were patterned by optical lithography. Volt- 
ages applied to a homogeneous top gate electrode were 
used to tune the electron density n e of the 2DEG. Trans- 
port experiments were carried out in a 3 He/ 4 He dilution 
refrigerator and in a standard 4 He cryostat at temper- 
atures between 60 mK and 20 K. The resistivity was 
measured in a standard 4-terminal setup, probing sample 
areas of 40/xm x lOO^im. For tlqd < 2 x 10 14 m~ 2 , the 
electron gases behave like ordinary low-mobility san*»les, 
in which the mobility is dominated by the SAQDs. til 
At higher dot densities of tiqd ~ 2.5 x 10 14 to~ 2 (i.e. 
the samples labelled by 14a, 14b and 24b, see inset in 
Fig. la), we find pronounced oscillations as a function 
of a magnetic field B applied perpendicular to the plane 
of the 2DEG. The oscillations are periodic in B, with 
a period of SB = 400 ± 20 mT in sample 14a, which 
increases to 560 ± 20 mT in sample 24b. The period 
is independent of n e (Fig. la). The relative oscillation 
amplitude Sp/p increases as n e is increased, with a max- 
imum of Sp/p ~ 0.12. Measurements in tilted magnetic 
fields show that the oscillations are determined solely 
by the magnetic field component perpendicular to the 
plane of the 2DEG, and have therefore a purely dynami- 
cal (no spin) origin. These oscillations are superimposed 
on a broad resistivity background, which vanishes around 
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B = 3T. We do not observe Shubnikov-de Haas oscilla- 
tions for filling factors larger than 4. 



23 
TEM 


13 


7M 


14a 


24b 


14b 


\ TEM 
\25 


15 




FIG. 



1. (a) 
2.37 x 



Magnetoresistivity 

io 15 ™- 2 ' 



sample 



14b for 



of 

■m 2 (D), 1.16 x 10 15 m- 2 («), and 
0.97 x 10 15 m -2 (A); the resistivity data for the two lower 
electron densities are scaled as indicated. Pronounced mag- 
neto-oscillations are observed for \B\ < 2T. Their period and 
phase is independent of n e . In addition, a positive magne- 
toresistivity is observed for \B\ < 0.5T, with the weak local- 
ization peak superimposed. The inset shows a sketch of the 
wafer under study, uqd increases from top to bottom (edge 
of the wafer), and has been estimated by TEM studies as 
n QD = 2.5 x 10 14 m~ 2 in sample 23, and n QD = 3.0 x 10 14 m~ 2 
in sample 25, respectively. The magneto-oscillations are ob- 
served on samples shown in dark gray, (b) Temperature de- 
pendence of the magnetoresistivity in samples 14b and 24b. 



In addition to a weak localization peak, we observe 
for large electron densities a well-pronounced minimum 
in the magnetoresistivity around B = 0, which extends 
up to a few hundred mT. In contrast to the magneto- 
oscillations, this minimum shows no significant temper- 
ature dependence up to about 2.5 K (Fig. lb). Figure 
2, shows additional observations that coincide with the 



presence of the features discussed above: (i) a fixed point 
is observed in the temperature dependence of p xx (n e ) 
at B = 0, separating a region of high electron density 
(n e > 1.2 x 10 15 m -2 ) with a metallic temperature de- 
pendence, from a low density—region showing an insulat- 
ing temperature dependence.113 (ii) the electron mobility 
fi shows a maximum as a function of n e (upper inset in 
Fig. 2). 
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FIG. 2. Metal-insulator transition as observed as a func- 
tion of n e and temperature in two samples. The upper inset 
shows the electronic mobility fi as a function of n e for two 
samples. Sample 24b shows a maximum in the electron mo- 
bility, in contrast to sample 24a. Lower inset: localization 
length vs correlation length j/o for the correlation function 
discussed in text in a ID system, and selected energy values. 
koc increases rapidly for larger yo, with a rate that depends 
on the energy E of the state. An even faster enhancement of 
lioc is anticipated for a 2D system. Algorithm to calculate li oc 
is described in Ref. 19. 

The presence of these effects depends on the cooldown 
cycle. Hence, the arrangement of charges in the po- 
tential landscape during cooldown is a crucial factor. 
However, we have been unable to establish a systematic 
dependence on the cooldown speed or on the biasing con- 
ditions during cooldown. The simultaneous occurrence 
of the transport features described here indicates a pos- 
sible unifying explanation which is developed below. 
In the following, we focus on sample 14a and interpret its 
behavior, starting with the periodic magneto-oscillations. 
The experiments clearly suggest a magnetic flux effect 
as its origin. It is well-known that Aharonov-Bohm 
oscillations with period AB = life A (A denotes the en- 
closed area) ensemble- average to zero, since the phase 
of each interference loop at B = is loop-specific.E-H la 
contrast, Altshulcr-Aronov-Spivak (AAS) oscillationslij 
always show zero phase shift at B = 0, and consequently 
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survive ensemble-averaging. We thus conclude that the 
observed oscillations are AAS oscillations with a charac- 
teristic enclosed area of 5.2 x 10~ 15 m 2 . 
Similar oscillations have been recently detected in hexag- 
onal antidot lattices, while no signature of such oscilla- 
tions has been found in rectangular lattices with other- 
wise comparable parameters Ji3 In Ref. [l6], the authors 
explain this observation in terms of semiclassical trajec- 
tories of the electrons around the antidots: in the hexago- 
nal array, the probability for the electrons to be scattered 
in loops around the antidots is greatly enhanced, leading 
to a significant AAS oscillation amplitude, in contrast to 
the more open geometry of rectangular lattices. 

Let us, based on this knowledge and our observations, 
assume that the potential landscape generated by the 
charged SAQDs and the strain fields reveals at least 
short range hexagonal ordering to a significant degree. 
The measured AAS period would then suggest an inter- 
dot separation of a « 78nm, corresponding to a density 
of SAQD points of 1.9 x 10 14 to~ 2 . This is in reasonable 
agreement with 2 x 10 14 to~ 2 < uqd £3x 10 14 to~ 2 , as 
determined in TEM studies, although iiqd as determined 
from resistivity studies is somewhat higher .EilO Further 
estimates are consistent with this picture. We have deter- 
mined the phase coherence length l^, from fits of the weak 
localization peak for p xx < h/e 2 and find 500nm 
at n e = 1.8 x 10 15 m~ 2 , and increasing slightly as n e 
increases. The electrons are thus clearly able to eclipse 
the area of the locally hexagonal ordered cells in a phase 
coherent fashion. As a further consequence, quantum 
point contacts are formed in between the lattice points, 
with a maximum of two occupied ID spin-degenerate 
subbands at the electron densities under investigation. 
This explains the quenching of the Shubnikov-de Haas 
oscillations above filling factor 4. Using these numbers 
and assuming that the hexagonal arrangement contains 
approximately circular scatterers of electrostatic radius 
tqd at the Fermi energy, we can estimate tqu w 20nm. 
Hence, 2tqd ~ f , which has been shown to be a partic- 
ularly effective geometry for -observing AAS oscillations 
in hexagonal antidot lattices.li3 We must emphasize that 
this analysis does not rely on truly long-range order, but 
rather only on the local hexagonal relative arrangement 
of the charged SAQDs. Nor do we know unambiguously 
that the ordering is really hexagonal: other types of 
correlations in the disorder potential, e.g. possibly a 
random trimer superlattice, may give rise to B-pcriodic 
coherent backscattering as well. Note, however, that the 
relative amplitude of the oscillations is comparable to 
that one observed in Ref. [ll| which suggests a similar 
degree of ordering. 

We proceed by constructing the relation between the 
AAS oscillations, the "metal- insulator transition" (MIT), 
and the maximum in fj,(n e ). If the density of SAQDs is 
low, there is no positional short-ranged order, and even 
homogeneous charging of the dots results in a potential 



landscape dominated by the disordered dot distribution 
of charged scattering centers. As the interdot separation 
is still large, the backscattering processes inherent to the 
AAS effect are unlikely, and no magneto-resistance os- 
cillations occur. Similarly, the mobility of the system 
is dominated by random scattering events from the dots, 
yielding an insulating behavior for all values of electronic 
density. As the density of dots increases, there is an ef- 
fective repulsive interdot interaction which facilitates the 
appearance of short-ranged Wider so as to relax the local 
elastic stress in the system.0il3 Upon formation of the 
2DEG, the dots so formed acquire charge with a distri- 
bution which is likely correlated to the loading state of 
nearby dots (notice that dot-dot interaction is quite large, 
as (2e) 2 /na w 6meV, for doubly-charged dots embedded 
in GaAs, k ~ 12). The short-range order yields enhance- 
ment of backscattering and AAS oscillations, as well as 
the appearance of relatively extended states due to the 
effectively weaker, more correlated, potential disorder. 
As the electron density increases, the AAS oscillations 
become more prominent, indicating increasing popula- 
tion of the corresponding states. The resulting electronic 
states are then capable of exhibiting a drastically differ- 
ent transport behavior as a function of electron density: 
insulating for low densities, and metallic for higher con- 
centrations. This change is the result of the Fermi level 
sweeping through the density of states, sampling ever 
more extended states. 

It is worth mentioning that our experiments show a sharp 
transition to such a correlated state as a function of uqd, 
given the fact that tiqd varies by 15 % between sample 
24a and 24b (or sample 13 and 14a) at most. Sample 
15 shows conventional behavior, which, however, may be 
due to its proximity to the wafer edge. 

The observed transport behavior is qualitatively con- 
sistent with the appearance of 'extended' states in the 
2DEG (with localization lengths much longer than a 
few SAQD spacings), and reminiscent of the situation 
recently analyzed in disordered but correlated ID sys- 
tems. Izrailev and coworkers have demonstrated that 
even in a ID system, potential correlations in a disordered 
structure can give rise to true mobility edges, where the 
localization length would diverge over a certain energy 
window.Ea Although the situation for 2D is not as firmly 
established, one expects that the higher connectivity of 
electronic paths in 2D would give rise to mobility edges 
even with weaker potential correlations than in a corre- 
sponding ID case^as numerical evidence in a variety of 
systems suggests.EJ 

In this context, it is instructive to estimate the po- 
tential correlations produced by the scattering centers 
(the charged dots) on the 2DEG carriers. The entire 
SAQD charged array generates a disordered Coulomb 
potential given by V(r) = J2j [F^~T e ~ 7 ^ r ~ Rj ^' wnere 
each dot has effectively trapped a charge qi and the bare 
dot potential is assumed screened by the 2DEG (the full 
potential surely includes also local strain components). 
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Correspondingly, the potential correlations are given by 
(V(x)V(x + y)), where the angular brackets denote an 
average over the dot population. For a random charge 
distribution, {qiqj) — Sijq 2 , and a random distribution 
of dots, the potential correlator has an asymptotic be- 
havior given by « Vge~ v ^ vo (a/y), where the correlation 
length j/o is a monotonic function of the screening length 
7 _1 , which in turn is a function of electron density in the 
system, and a is the mean dot separation. One expects 
that as the electron density increases and provides bet- 
ter screening, the characteristic yo also increases, result- 
ing in stronger (longer range) correlations in the system. 
As this expression neglects the built-in effective interac- 
tion between dots introduced by the strain fields, as well 
as the correlations in potential strength of the charged- 
dot landscape, this result provides a pessimistic estimate 
of the true correlator in our system. Nevertheless, even 
such a 'weak' exponential correlator in a ID system can 
be shown to result in increasing localization lengths fpc 
larger yo values, as shown in the lower inset of Fig. 2x3 
The rather strong potential correlations that exist in our 
2D system are then anticipated to yield long localization 
lengths which clearly exceed the phase coherence length 
of the system, and explain the complex behavior we ob- 
serve. 

Within this interpretation, the character of the MIT ob- 
served in our system is quite different from those reported 
in other systems. One striking feature of the MIT ob- 
served in, e.g., Si-MOSFETs is its suppression in mag- 
netic fields applied parallel to the plane of the electron 
gas. This suppression has been attributed to complete 
spin alignmentO This is in sharp contrast to our system, 
where the MIT is present up to the strongest magnetic 
field accessible to us, i.e. 12 T. 

A positive magnetoresistivity around B = is also 
present in the metallic regime. This resistivity minimum 
is robust as the temperature is increased (Fig. lb), in- 
dicating a classical origin. Furthermore, the minimum 
vanishes as n e is reduced into the insulating regime (Fig. 
la). Similar effects are well-known from electron gases 
with double subband occupancy:E3 due to the different 
densities and mobilities of the electrons in the two sub- 
bands, a positive magnetoresistance evolves. Since we 
cannot determine the individual electron densities in the 
two types of states, we cannot fit p X x{B). Qualitatively, 
however, it is known that this effect is strong if the differ- 
ence in the mobilities is significant, and consistent here 
with the presence of a 'band' of high mobility (extended 
states) next to another 'band' with low mobilityHa 

In conclusion, the following qualitative picture 
emerges: the presence of AAS oscillations indicates some 
degree of short-range ordering, possibly in a hexagonal 
geometry, of the potential landscape induced by the InAs 
SAQDs. This leads to the formation of a different class 
of states, which may be seen as an extended-state "mini- 
band" arising from the local ordering in the dot potential. 
These states become populated as n e is increased above 
a threshold density, which corresponds approximately to 



the fixed point of the MIT. For higher electron densi- 
ties, scattering between two bands with quite different 
mobilities can occur, which reduces the overall mobility. 
The positive magnetoresistivity around B — indicates 
such a scenario. Our data thus represent the first trans- 
port signature that two-dimensional, correlated potential 
landscapes can be formed by self-ordering of charge be- 
tween self- assembled quantum dots, and raises several 
questions about mobility edges and miniband design by 
self-ordering in two dimensions, which will hopefully be 
addressed in future experiments and theoretical studies. 
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